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_i ■ Several dark energy models on the brane are investigated. They are compared with corresponding 

theories in the frame of 4d Friedmann-Robertson-Walker cosmology. To constrain the parameters 
of the models considered, recent observational data, including SNIa apparent magnitude measure- 
ments, baryon acoustic oscillation results, Hubble parameter evolution data and matter density 
perturbations are used. Explicit formulas of the so-called state-finder parameters in teleparallel 
theories are obtained that could be useful to test these models and to establish a link between 
Loop Quantum Cosmology and Brane Cosmology. It is concluded that a joint analysis as the one 
developed here allows to estimate, in a very convenient way, possible deviation of the real universe 
a* cosmology from the standard Friedmann-Robertson-Walker one. 

I. INTRODUCTION 

^ 1 A number of difficult problems in cosmology have been put forward by the discovery of the accelerated expansion of 
the universe 0, H[ . This cosmic acceleration can be explained via the introduction of dark energy (for a recent review, 
see [3] ) . It follows from recent observational results that dark energy currently accounts for about 73% of the total 
mass/energy of the universe [j. It appears to have rather strange properties, as a negative pressure and/or negative 
entropy, the fact that it is undetectable in the early universe, and so on. It is still not excluded, however, that General 
Relativity (GR) and the ensuing vacuum fluctuations (as those leading, e.g., to the Casimir effect) could lead to an 
explanation of the issue (see e.g. [1, §)■ One should also stress the following important connection: with the help of 
an ideal fluid, GR can actually be rewritten, in an equivalent way, as some modified gravity (for a recent review, see 

0). 

For dark energy with density pd and pressure pu, the equation of state (EoS) parameter wu, 



(N 



x 



w-d = Pb/pb < , (1) 

is known to be negative. Moreover, astrophysical observations clearly favor, up to now, the standard ACDM cosmology. 
Dark energy as just a cosmological constant (wd = —1) is the simplest and maybe best preferred model from the 
theoretical viewpoint, too. In this model over 70% of the current energy budget is dark energy (Einstein's cosmological 
constant A), in perfect accordance with the data coming from observations, as reported above. 

Presently several independent observational procedures provide strong evidence in favor of the ACDM model, in 
particular SNIa apparent magnitude measurements as a function of the redshift, cosmic microwave background (CMB) 
anisotropics, baryon acoustic oscillation (BAO) peak length scale measurements, Hubble parameter determinations 
etc., but the error bars associated with all these classes of data are still too large to allow for a significant observational 
discrimination between the ACDM model and other existing, theoretically well grounded alternatives to it. 
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When w < — 1 (phantom dark energy) [7| , we face up the most interesting and less understood theoretical possibility. 
A simultaneous violation of all four energy conditions occurs in this case and the involved field is unstable, although 
it could perhaps be made stable in classical cosmology [8(. For a universe filled with phantom energy there are 
many possible new scenarios for the end of such universe, among which the most typical cases are those of a Big 
Rip singularity fA l9l-[l2j and of a sudden future singularity [l3j . However, a final evolution without singularities is 
also possible: if the parameter w asymptotically tends to —1, and the energy density increases with time or remains 
constant, no finite-time future singularity will be ever formed 

EMI, 

also, if the universe starts to decelerate in 
the far future. In any case, if the energy density grows up to some threshold value, the disintegration of any bound 
structure will eventually occur, in a way quite similar to the case of the Big Rip singularity, but this may only happen 
faraway in the future evolution. The dark energy pressure is expressed as a function of the density, as 

p = g(p)- (2) 

In this paper a number of dark energy models on the brane will be considered (for a general introduction to brane- 
world cosmology, see The theoretical predictions of these models will be compared with the various types of 

existing independent data observations, including the luminosity distance modulus vs redshift for SNe la, the data 
accumulated on the evolution of the Hubble parameter H(z), the latest baryon acoustic oscillation (BAO) results, and 
matter density perturbation data. In Sect. II a brief overview of the EoS fluid formalism is presented. A comparison of 
dark energy in Friedmann- Robertson- Walker (FRW) cosmology and on the brane is carried out. The main constraints 
coming from the observational survey data will be analyzed in Sect. III. In the following two sections, Sects. IV and 
V, we will study the simplest ACDM model on the brane and will show that a careful joint analysis of the various 
observational data allows to estimate, in a clear fashion, any possible deviation of our cosmological model from the 
standard FRW cosmology. Due to the increasing interest in teleparallel theories (F(T) models), we deal with them 
in Sect. VI where we find explicit formulas for the so-called state-finder parameters in such theories, which could be 
useful to test the models proposed there, in particular to test Loop Quantum Cosmology and Brane Cosmology, and 
specific relations between the two. Finally, Sect. VII is devoted to conclusions. 



II. COMPARISON OF DARK ENERGY IN FRW- AND BRANE-COSMOLOGY 



We start with a brief description of dark energy models in the frame of the FRW cosmology. The cosmological 
equations corresponding to a spatially flat universe, endowed with a metric 

ds 2 = -dt 2 + a 2 {t)(dx 2 + dy 2 + dz 2 ), (3) 

are the following 

:)'-§■ >->(!>+»>■ w 

with p and p, respectively, the total energy-density and pressure, while a is the scale factor, the dot means time 
derivative, and natural system units are being used, with 8nG = c = 1. 
For dark energy, the EoS can be rewritten, for convenience, in the form 

Pd = -Pd - /(pd), (5) 

being /(pd) a function of the energy-density. We observe that /{pd) > corresponds to w < — 1, and J{pd) < to 
w > —1. The future evolution of the universe depends on the EoS for dark energy chosen. Let us here describe two 
main cases. 

(i) Evolution without future singularities. This case includes a so-called "Little Rip" (these models are described in 
detail in [Til UH). The dark energy density grows with time so slowly that a Big Rip cannot occur in finite time. For 
the realization of this scenario, one needs the asymptotic behavior of the function to be g{pv) ~ Px>, as /? < 1/2. 
Eventually, a dissolution of all bound structures will also take place in the future. 

One should note that, for some specific equations of state with branch points, a de(phantomization) process can 
occur 20]. Therefore, after the acceleration phase, a slowdown of the future universe might be possible. In other 
words, the universe may be decelerating in the future. 

Interesting alternatives to the ACDM model are models in which the dark energy density tends asymptotically to 
a constant value (an "effective cosmological constant" [Hj])- One should remark that if the value of this "effective 
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cosmological constant" is sufficiently large (for example, if A e // ~ 1 in Planck units) a possibility of disappearance of 
bound structures due to the enormous acceleration of the universe still remains. 

(ii) Evolution with finite-time singularities. If g(po) ~ Pp, with (3 > 1/2, the dark energy density grows so rapidly 
that the universe ends its existence in a singularity of Big Rip type or in a type III singularity, according to the 
classification in pl| . The key difference between these singularities is that the energy density in the second case grows 
so rapidly with time that the scale factor does never reach the infinite value. These scenarios can be realized only in 
the case of having a phantom energy. Another interesting case occurs if /(pd) — > ±oo at pn = PDf, i-e., the dark 
energy pressure becomes infinite at finite energy density. The second derivative of the scale factor diverges, while the 
first derivative remains finite. 

As an alternative to the FRW cosmology let us consider the simplest brane model in which spacetime is homogeneous 
and isotropic along three spatial dimensions, being our 4-dimensional universe an infinitesimally thin wall, with 
constant spatial curvature, embedded in a 5-dimensional spacetime [H, [23[. In the Gaussian normal coordinate 
system, for the brane which is located at y — 0, one has 

ds 2 = -n 2 dt 2 + a 2 (t, y)^ ij dx i dx j + edy 2 , (6) 

where 7^ is the maximally 3-dimensional metric, and e = 1, if the extra dimension is space-like, while e = — 1 , if it is 
time-like. Let t be the proper time on the brane (y = 0), then n(t, 0) = 1. Therefore, one gets the FRW metric on 
the brane 

dsf y=0 = -dt 2 + a 2 (t, 0)j l3 dx l dx 3 . (7) 
The 5-dimensional Einstein equations have the form 

Rab - t;9abR = X 2 Tab + ^9ab, (8) 

where A is the bulk cosmological constant. The next step is to write the total energy momentum tensor Tab on the 
brane as 

T£ = Sp(y), (9) 

with Sg = di&g(—pi,,Pb,Pb,Pb, 0), where pb and pb are the total brane energy density and pressure, respectively. 

One can now calculate the components of the 5-dimensional Einstein tensor which solves Einstein's equations. One 
of the crucial issues here is to use appropriate junction conditions near y — 0. These reduce to the following two 
relations: 

ndy { 0+ 3 2 ady\ 0+ 6 



A k C 



After some calculations, one obtains the following result 

H 2 = e X ^ + y-- + - 4 . (11) 
3d d a a 1 

This expression is valid on the brane only. Here H = d(t,0)/a(t,0) and C is an arbitrary integration constant. The 
energy conservation equation is correct, too, 

Pb+3% b +p b ) =0. (12) 
a 

Now, let pb = p + A, where A is the brane tension. For a fine-tuned brane with A = eA 2 x 4 /6, we have the equation 
(for k = 0) 

In what follows we will consider a single brane model which mimics GR at present but differs from it at late times. 
We set 87rG = ecrx 4 /6. As can be checked, two cases arise: e = 1 and A > 0, and e = — 1 and A < 0. For simplicity, 
we also set C — (the term with C is usually called "dark radiation"). In fact, allowing C 7^ does not lead to 
additional solutions on a radically new basis, in the framework of our approach. Eq. (1131) can be simplified to 

a 2 p A. , p 
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Eq. (fl~4|) . for p « [AL differs insignificantly from the FRW equation. The brane model with a positive tension has 
been discussed in |24l426l | in the context of the unification of early- and late-time acceleration eras. The braneworld 
model with a negative tension and a time-like extra dimension can be regarded as being dual to the Randall-Sundrum 
model [27M29I which, as we will show in Sect. VI A, is in fact equivalent to the effective formulation of loop quantum 
cosmology. Note also that, for this model, the Big Bang and Rip singularities are absent (see Sect. VI A). And this 
fact does not depend upon whether or not matter violates the energy conditions [30l ]. This same scenario has also 
been used to construct cyclic models for the universe [3lj |. 

One can assume that in our epoch p/2X << 1 and thus there is no significant difference between the brane model 
and FRW cosmology. But the universe evolution in the future, for brane cosmology, can in fact differ from such 
convenient cosmology, due to the non-linear dependence of the expansion rate on the energy density. 

The EoS formalism for dark energy models on the brane was considered in [32j]. Here we briefly describe this 
approach. One gets the following link between time and dark energy density, assuming that pD » p rn : 



1 f PD do 

t(p D ) = -— ^ • (15) 

For the present time, to, we can set to = 0. For the scale factor as a function of the dark energy density, we have the 
same relation as in the FRW cosmology, namely 

fl f PD dp \ 

a = a Q exp - / — — . (16) 



, 3 Jpno 9{p) 

In the case of a positive tension, the following possibilities can be realized: 

1. If the integral (fT5")) converges while (|16p diverges, we have a Big Rip. It is interesting to note that the Big Rip 
on a brane considered in [33| occurs faster than in the FRW cosmology. 

For the simplest EoS with constant state parameter wq = —1 — a 2 , the function g{po) — o?pd- If Pd » A, 
then the dark energy density grows with time substantially faster than in ordinary cosmology (A — > oo). 

2. If the integrals (|15p and ([To]) diverge when po — > oo, then a Little Rip occurs. The acceleration of the universe 
increases with time definitely faster than in the FRW universe, owing to the brane tension (see the corresponding 
time equation for the case (1)). 

3. Asymptotic de Sitter expansion is realized if g — ¥ for pd — > pDf, and the integral (|15[) diverges. 

4. There is a type III singularity if both integrals converge when pjj — > oo. 

5. If g(pD) °° for Pd — ► PDf, the universe ends its existence in a sudden future singularity. 
The case of negative tension allows for the following interesting possibilities: 

1. Asymptotic de Sitter expansion, if g(po) for pD — > PDf- 

2. An asymptotic breakdown (i.e. the rate of expansion of universe tends to 0) will occur if g{po) — > for pD — > 2A. 

3. A sudden future singularity, if /(pd) oo when pd — > PDf- 

One should note that dark energy with an EoS such that f(p) ~ p 1 , with 7 < 2, leads to a Big Rip on the brane 
while, in the case of the conventional FRW universe, such dark energy leads to a Little Rip only. 



III. OBSERVATIONAL DATA 



The parameters of the cosmological models can be determined from a strict comparison of their predictions with 
accurate observational data. We here consider the data coming from SNe observations, the evolution of the Hubble 
parameter, baryon acoustic oscillation, and the evolution of matter perturbations. 
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A. SNe observations. 



The modulus p vs redshift z = a$/a—l relation corresponding to type la supernovae from the Supernova Cosmology 
Project [Hj], is, as well known, 

p(z) = pa + 51gL> L (z). (17) 
The relation for the luminosity distance Dl(z) as a function of the redshift in the FRW cosmology (FC) is 

D[ c = ^-(l + z) [ Z h- 1 (z)dz, h(z)=[n m o(l + z) 3 + Q DO F(z)} 1/2 (18) 
■"o Jo 

Here, f2 m o is the total fraction of matter density, Qdo the fraction of dark energy energy density, and Hq the current 
Hubble parameter. The constant value po depends on the chosen Hubble parameter: 

p = 42.384 - 5 log ft, h = i7 /100km/s/Mpc 

The function F{z) — pd{z) / Pdo can be determined from the continuity equation 

PD -3-g(p D )=0, (19) 
a 

which can be rewritten as 



g(y) 



-31n(l + z). (20) 



For simplicity, we neglect the contribution of radiation. 

For cosmology on the brane (BC), Eq. (|T5)) can be rewritten as 

Df c = -^(1 + z) f h-\z) [1 + Sh 2 (z)] - 1/2 (1 + 5) l ' 2 dz (21) 
H a Jo 

where for convenience the parameter S = po/2X has been introduced. For the analysis of the SNe data one needs to 
calculate the parameter x 2 , which is defined by 

2 _ \Pobs(Zj) ~ PthjZj)} 2 . . 

xsn — 72 ' y zz ) 

i * 

where Oi is the corresponding la error. The parameter po is independent of the data points and, therefore, one has 
to perform a uniform marginalization over pq. Minimization with respect to po can be done by simply expanding the 
Xsn with respect to po, 

X 2 sn = A- 2p B + p 2 C, (23) 

where 

[pobs(zi) - pth(zi;po = 0)] 2 



B _ Pobs{Zj) - Pth(Zj) _ J_ 

The expression f|23p has a minimum for /io = B/C at 

Xsn = A - B 2 /C. 

One can minimize Xsn instead of Xsn- Following (36|, one determines the 68.3% confidence level (CL) by A% 2 = 
X 2 ~ Xmin < 1-0 f° r the one-parameter or 2.3 for the two-parameter model. Similarly, the 95.4% confidence level is 
determined by A\ 2 = X 2 ~ X 2 m in < 4.0 or 6.17 for the one- and two-parameter models, respectively. 
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z 


H obs (z) 






km s _1 Mpc -1 


km s^ 1 Mpc -1 


0.090 


69 


12 


0.170 


83 


8 


0.270 


77 


14 


0.400 


95 


17 


0.480 


97 


62 


0.880 


90 


40 


0.900 


117 


23 


1.300 


168 


17 


1.430 


177 


18 


1.530 


140 


14 


1.750 


202 


40 



Table I: Hubble parameter versus redshift data from [371 ]. 
B. Hubble parameter. 

The evolution of the Hubble parameter with time in the past is now well observed. The Hubble parameter depends 
on the differential age of the universe as a function of the redshift, in the form 



dt = 



1 dz 
Hl + z' 



Therefore, a determination of dz/dt directly measures H(z). This is made possible through data we have on the 
absolute age of passively evolving galaxies, determined from fitting stellar population models. We use the 11 datapoints 
for H(z) from [37| for constraining the model parameters. These data are listed in Table I. The theoretical dependence 
of the Hubble parameter in the brane model is 



H{z) = H h{z){l + Sh\z))^ 2 {l + S)- 1 



-1/2 



The parameter \h ^ s here 



Xh 



£ 



[H bs(zi) - H th (zi)Y 



(24) 



(25) 



We need to perform a uniform marginalization over the parameter Hq. Again, we can expand 

x \ =A 1 -2BiH + H%C 1 , 



H bs{ziY 



E(zi)H obs (zi) 



The parameter x% has a minimum at the point Hq = B\jC\, 

xl = A 1 -Bl/C 1 . 
As in the case of the SNe data, we could minimize Xh instead of Xjj ■ 



C. BAO data. 

To constrain cosmological parameters using BAO data we follow the procedure described in [38j]. We use the 
measurements of the acoustic parameter A(z) from [38(, where the theoretically predicted A t h(z) is given by the 
relation 



D v (z)H y/Q m0 

Ath(z) = , (26) 
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where Dy{z) is a distance parameter defined as 



D v {z) 



Here, cIa(z) is the angular diameter distance 



(l + zfd\(z)- C ' 



H{z) 



1/3 



(27) 



d A (z)= V \ Z) y( z )=f Z ^- E(z) = H(z)/H . (28) 

H (l + z) Jo £(z) 



Using Eqs. ([26 ]) -([28 ]) , we have 



= V^o 4w~v (29) 
and using the WiggleZ ^4 & s (z) data from Table 3 of [38j |. we compute Xa to be 

xi - AA T (C A )- 1 AA. (30) 

Here, AA is a vector consisting of differences, AAi = A t h(zi) — A b s (zi) and C^ 1 is the inverse of the 3x3 covariance 
matrix given in Table 3 of [38| . 

D. Matter density perturbations. 

As was shown in [3^] one can neglect the density perturbations of dark energy. In this case the dark matter 
perturbations effectively decouple from DE perturbations. The equation that determines the evolution of the density 
contrast S in a flat background filled by matter with density p m is 

S m + 2HS m = -p m S m - (31) 

It is convenient to introduce the growth rate function of the perturbations / = c?ln(5 m /dlna. Using the FRW 
equations, one gets the following equation for / 

^ + / 2 + (^ + 2)/-^ m (l + P/2A)-^0, (32) 

where Cl m is the matter fraction of the total energy-density fi m = O m o(l + z) 3 /[^mo(l + z) 3 + ^IdqF(z)]. Finally, 
using the relation 

d ^ d 



d In a dz 



and taking into account that 



we get 



a a 2 p + P 



(1 + P/A), 



- (1 + z)f z + f + (2 - + / - ^(1 + P/2A)- 1 = 0. (33) 

where fljj = po/ P and we have introduced the parameter A = (1 + p/A)(l + p/2A)~ 1 . 

For a dark fluid with given EoS, one can find the DE density as a function of the redshift z. Then, Eq. (|3"3"|) can be 
solved numerically. The observational data for the growth factor f b s at various redshifts are given in Table II. 
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z 


fobs 


Ref. 


0.15 


0.51 ±0.11 


[40], [41] 


0.32 


0.654 ±0.18 


[42] 


0.35 


0.70 ±0.18 


[43] 


0.55 


0.75 ±0.18 


[44] 


0.77 


0.91 ±0.36 


[45] 


1.4 


0.90 ±0.24 


[46] 


3.0 


1.46 ±0.29 


[4JJ 



Table II: Available data for the growth factor at various redshifts from the change of the power spectrum Ly-a forest data 
in SDSS. 



Data sets 


SNe 


H 


BAO 


SNe±H±BAO 


SNe±H±BAO±F 


Xmm 


5 = 


„ 7 „ o -0. 019, -0.039 
lzz +0. 020, +0.039 


n _„_-(!. 089. -0.213 
U. IOO +0 .066.+0 .117 


n CQQ-0. 028, -0.059 
U.Oyy +0 .026, +0.049 


„ 71 „-0. 011, -0.027 
U -' iZ +0. 016, +0.031 


„ 71 9 -0. 009, -0.025 
u -' iz +0. 018, +0.033 


562.39 


8 = 0.05 


„ .-0.016,-0.036 
' * 4 +0. 016, +0.034 


„ 7s7 -0. 058, -0.130 
' 8 '+0. 047, +0.086 


fl oqc— 0.031,— 0.065 
U.D5D +0 .028, +0.054 


fl 7OQ-°' 012 > -0.028 
U - ' Zy +0. 013, +0.028 


fl 7OQ-0. 012, -0.027 
U - ' Zy +0. 013, +0.027 


565.91 


8 = 0.10 


„ 7 r Q -0. 013, -0.030 
U. /oy +0 .014.+0 .032 


fl »1 q-0 046, -0.102 
"■oJ-J+0. 039, +0.072 


fl C7O-0. 028, -0.072 
U-D'^+0. 031, +0.059 


fl 7/|q-0.014,-0.028 
U. dO +0 .on,+o .025 


fl 7/|q-0.014,-0.028 
U. dO +0 ,oio,+0 .023 


571.38 



Table III: Best fitting values for Oa within la and 2a errors for various 8 from observational data analysis for the ACDM model 
on the brane. In the last column the minimal value of Xmin — Xsjv ± Xlf ± Xa ± X/ is given. 



IV. ACDM MODEL ON THE BRANE 



First, we consider the very simple cosmological model on the brane with vacuum energy po = A = const. This 
model coincides in the future with the FRLW cosmology with a redefined cosmological constant. The asymptotic 
behavior of the scale factor is 

a(t) ~ ao exp (A^/3) 1 / 2 * L t^oo, A eff =A(1 + A/2A). (34) 

One can consider the ACDM model on the brane as a one-parametric model, at fixed values of 5. The results of 
the calculations corresponding to this case are given in Table III (we have also included in our considerations the case 
5 = 0, i.e. FRW cosmology, for comparison). The BAO data favor smaller values of Udo than the H(z) and SNe 
data. The optimal value of 51a (that is, ft d when po = A) is closer to the one coming from the SNe data analysis 
only. One easily sees that the addition of the observational data for the matter density perturbations does not change 
the best-fit value of fl\ from the SNe±H±BAO analysis. One can also conclude that the best consistent description 
of all observational data is realized in the frame of the FRW cosmology (<5 = or A — > oo). The minimal value of the 
total x 2 is 562.39. As 5 grows the corresponding \ 2 increases. 

For 5 > we have the following picture. The analysis of the data sets does not yield separately a significant 
constraint on the maximal value of 5. The parameter x 2 for SNe, BAO and matter density perturbation data grows 
very slowly with increasing 5. For instance, for 5 = 0, XsNmin = 553.18, while for 5 = 0.1 the minimal value of Xsn 
is 553.34. The data on the evolution of the Hubble parameter are more sensitive to increasing 5: for 5 = 0.1 we found 
that Xh min — 8.12, in comparison with x% m in = 7-62 in the FRW model. But one can see that, at 5 ~ 0.05, the 
la intervals of the possible values of il\ for SNe and BAO data do not intersect. The 2a intervals for from these 
data sets do not have common points for <5 w 0.10. Therefore, one can estimate the maximal value of 5 from the joint 
analysis of all observational data sets. 

For doing this, we consider the ACDM model on the brane as a two-parametric one, with free parameters 5 and 51a- 
One can see that, although the areas corresponding to the la and 2a contours from the SNe, BAO and H(z) data 
analysis are sufficiently large (see Fig. 1), these contours intersect in a quite narrow region of the parameter space. 
Joint data analysis allows us to define the la and 2a contours in the J7a - 8 parameter space (Fig. 2). Therefore, we 
can estimate the upper limit of the parameter 5 at which the ACDM model is relevant to the observational data. 

V. OTHER DARK ENERGY MODELS ON THE BRANE 

Let us now consider the following model with a quite simple EoS, 



i(Pd) = a 2 p Da ( — ) , (35) 

\PD0j 
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0.04 - 



0.02 



Figure 1: The 68.3% (upper panel) and 95.4% (bottom panel) confidence level contours in the Q,a 
the analysis of SNe (solid line), H(z) (bold solid line) and BAO (dotted line) data, respectively. 



5 parameter space from 



where a and /3 are dimensionless constants. If j3 = 1, this reduces to an ordinary phantom energy model with 
constant EoS parameter w = — 1 — a 2 . From Eqs. (|15[) and (ITBl one can see that, for various values of /3, the model 
(|35f describes three types of future universe evolution: 

(a) Little Rip, for (3 < 0, 

(b) Big Rip, for < ft < 1, and 

(c) a type III singularity, for (3 > 1. 

Note that in the FRW cosmology the Little Rip occurs for f3 < 1/2. Simple calculations allow us to obtain the 
function F(z): 



F(z) 



f [1 - 3a 2 



;i-/3)ln(l + z)]^ , 
(l + z)- 3 " 2 , (3 = 1. 



(36) 



We consider the case f3 = (Little Rip) and 0=1. 

The analysis of observational data for the Little Rip model leads to the same conclusions as in the case of the ACDM 
model: when the brane tension decreases, the common area of the confidence level contours for SNe, H(z), and BAO 
data in the a 2 - floo parameter space decreases too, that is, the agreement with observational data becomes worse. 
In Fig. 3 the la confidence level contours from the data set analysis are shown. The results for the joint observational 



10 



6 



u.uo 


1 1 1 1 1 1 1 1 




07 




■ \ 




OR 






05 






04 






U.UO 






0.02 








0.01 











1 1 / 1 1 1 1 / 


1 1 ■ 



0.67 0.68 0.69 0.7 0.71 0.72 0.73 0.74 0.75 0.76 



Figure 2: The 68.3% (solid) and 95.4% (dotted) confidence level contours in the 51a - 5 parameter space from the 
SNc+ff (z)+BAO+matter density perturbation data analysis. The best-fit parameters for the observational data are ^a = 
0.712, 5 = with x 2 = 562.39. 



data analysis are depicted in Fig. 4. One can see that, for large values of S, the description of the observational data 
is better for larger a 2 . 

Similar results can be derived for the simplest phantom model with constant EoS parameter w . The best- fit 
parameters for the FRW cosmology are (w = —1.01, £Ido = 0.713) with \ 2 — 562.16 (recall, for comparison, that in 
the ACDM model we have the slightly larger value of Xmin — 562.38). In Fig. 5 the results of the joint analysis are 
depicted. For 5 > the observational data analysis speaks evidently in favor of w < — 1. 
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Figure 3: The 68.3% confidence level contours in the Qa - a 2 parameter space from the analysis of SNe (solid), H(z) (bold 
solid), and BAO (dotted) data for the FRW cosmology (A), 8 = 0.05 (B), and 8 = 0.10 (C), in the case of the Little Rip model 
with (3 = 0. We see that, for each of the data sets, the best-fit parameters correspond to the ACDM model (a 2 — 0). 
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Figure 4: The 68.3% (solid) and 95.4% (dotted) confidence level contours in the - a 2 parameter space from the analysis 
of SNe+//(,z)+BAO+matter density perturbation data for the FRW cosmology (A), 8 = 0.05 (B), and 8 = 0.10 (C), in the 
case of the Little Rip model with /3 = 0. The best-fit parameters are (a 2 = 0.03, Q D o = 0.713) with X 2 = 562.20 for the FRW 
cosmology, {a 2 = 0.08, Q D0 = 0.714) with X 2 = 563.73 for 8 = 0.05, and {a 2 = 0.10, Q D0 = 0.722) with X 2 = 566.78 for 
S = 0.10. 
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Figure 5: The 68.3% (solid) and 95.4% (dotted) confidence level contours in the Qa - wo parameter space from the analysis 
of SNe+//(z)+BAO+matter density perturbation data for 8 — 0.05 in the case of the simplest phantom model. The best-fit 
parameters are (wq = —1.08, Qdo = 0.718) with \ 2 = 563.84. 



VI. STATE-FINDER PARAMETERS IN TELEPARALLEL THEORIES 



In this Section we will compute the so-called state-finder parameters in universes described by teleparallel models 
which may be considered as an example of Loop Quantum Cosmology. The first of these parameters is the effective 
u! parameter and the second the deceleration parameter. Both are well-known in the literature. The other two 
parameters, introduced due to increase in the accuracy of cosmological data, were given in [48[ (see also |49j). with 
the purpose to advance beyond the effective w and deceleration parameters. The four parameters are defined as 
follows: 



1. The effective uj parameter is 



2. the deceleration parameter 



3. the jerk parameter 



4. and the snark parameter 



a 



J e// = -i-|^; (37) 



(38) 




H 3 H 3 H 2 



J -I 



3(q dec - 1/2) ' 



(40) 



We consider again a universe filled by a perfect fluid with EoS p = g(p) = — p — Hp)- Teleparallel theories in flat 
FRW cosmology are defined via a Lagrangian of the form Ct — VF(T) — Vp (see |50l- [52f for a review of the topic), 
where V is the volume of the spatial part and T — ~6H 2 is the so-called scalar torsion^m, 1531] . 

From this Lagrangian we can see that the conjugate momentum of V is given by py — j& = —4HF'(T), and thus 
the Hamiltonian is 

H = Vpv -£ = [2TF'{T) - F(T) + p] V. (41) 

It is well-known that in general relativity the Hamiltonian is constrained to be zero. This constrain leads to the 
modified Friedmann equation 

p=-2F'{T)T + F{T) = G{T), (42) 
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which is a curve on the plane (H, p) . 

Conversely, given a curve of the form p = G(T) for some function G, it could be obtained from the modified 
Friedmann equation, by choosing [52[ 

F(T) = -^f J §j=dT. (43) 

The modified Raychaudhuri equation is obtained from the modified Friedmann equation by taking its derivative with 
respect to time and using the conservation equation p — 3Hf(p), from where H = —jfpj^- Then, the dynamics of 
the universe is given by the modified Raychaudhuri equation and the conservation equation, i.e. by the system 

Hp) 

4G'(T)' (-44) 
3Hf(p), 

provided the universe moves along the curve p — G(T). 

To compute the state-finder parameters, one has to use the modified Friedmann, Raychauduri and conservation 
equations, to get 

w m _ , W(T)) ( 3f(G(T)) 
w e//(Tj — -1 jx^JT)' qdeciT) -~\2TG'(T) +1 

m = ^[^3 [(f(G(T)) + 1)(G'(T)) 2 - G"(T)f(G(T))] + 1, 

f(G(T)) [(f'(G(T)) + l)(G'{T)f G"(T)f(G(T))} 
1 ' (G'(T)) 2 f{G(T))+TG'(T) ' 1 ; 

These formulas mean that the parameters are functions of H. Since H (the current value of the Hubble parameter) is 
well-known, one can test all the F(T) models with current observations without actually solving them. In particular, 
as we will see, one can test loop quantum cosmology or brane cosmology in the Randall-Sundrum scenario, i.e., branc 
cosmology described in Sect. II. In fact, Hq could be calculated from measurements of the luminosity distance Dl{z), 
via the well-known formula (|48j|) 



H(z) 



d (D L {z) 



(46) 



dz \ 1 + z 

being z the redshift, or either as we have already explained in Sect. Ill B. 

A remark is important: the formulas (|45[) could be expressed in function of p via the relation T = G^ 1 (p). This 
will be relevant when we deal with loop quantum cosmology and brane cosmology. These expressions read 

f(p)(G~ 1 )'(p) ( 3f(p)(G-i)'(p) 

"e//(p)=-l- G . X{p) , Qdec(p) = - y 2G~Hj)) 

i{p) = ^r^y [(f'(p) + iKG-'Yip) + f(p){G~ 1 )"(p)] 4 



Hp) 

f( P )(G-iy(p) + G-i( P y 



*(P) = M , ^tJ if'iP) + ^(G-'Yip) + MiC- 1 )"^)}. (47) 



We start calculating the parameters for the simplest but one of the most interesting EoS, when the dependence 
between pressure and energy density is linear, i.e., for 

p = uip <^ f(p) = -(1 + ui)p. 

In that case, one has 

We// (T) = -l + (l + w)^ 7 , q dec (T) = - ( -hi ( ' 



TG' " v ' V 2 TG 
3{T) = - |(1 + ^T^f (-w(G') 2 + (1 + u)GG 



._. fl + u\G\-uG' + (1 + lj)GG"] . . 
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As an application, we study FRW cosmology with a small cosmological constant A. In this case G(T) 
from (|4"5|) . a simple calculation yields 



-j — A and, 



^e//(T) = W + 2(1 + W)-, grfec(T) 



3, 



j(T) = -(l+w)« 



2A 



+ 1, a(T) = u 1- 



2A 



1 



(49) 



We should remark that, when cj > —1, at late times p — > and thus T — > — 2A. Obviously one has w e //(T) — » — 1, 
1dec{T) — > — 1, — > 1 and s(T) — > 0. On the other hand, when w < — 1 at late times one has a Big Rip singularity. 

Recently, another interesting model in FRW cosmology has been introduced in [54| in order to deal with non-singular 
universes. One considers once again the curve p = G(T) = — -j — A, but with a non-linear EoS 



p{p) = - P -^f{p) = -p(l-P- 



Pi 



Pi 



(50) 



where pi is a constant satisfying A <C pi 



This model has two de Sitter solutions Hf = \ / 4 and H, 



, ...... ...... — . ., , and shows a universe evolving from an early 

inflationary phase (the de Sitter phase Hi) to late time accelerated expansion (de Sitter phase Hf), going trough 
a matter dominated phase which allows the formation of structures. It could be also viewed as a universe with a 
huge cosmological constant pi at early times, which evolves, at late times, towards a small cosmological constant A 
responsible of the current cosmological acceleration. 

In this case, for the current (small) value of T, using (|4"T|) . (|50l) and G~ 1 (p) = ~2(p + A), one gets 



eff(p) = -1 



J(P) 



p{Pi - P) 



Pi(p + A) 
9p 2 (p t -p) 

pKp + a ) 



qdec(p) = -1 - 



3p(p» - p) 



S{P) 



2pi{p- 
P 2 {P 1 - 



-A)' 
P) 



(51) 



A. Loop Quantum Cosmology and Brane Cosmology with a small cosmological constant A 

In loop quantum cosmology, which reminds of brane cosmology, the effective Friedmann equation depicts the 
following ellipse (see, for example, [521 ]) 

on the plane {H,p), where p c is the so-called critical density, that satisfies A <C p c - This curve can be written in two 
pieces, p m — G-(T) and p m = G+(T), where 



feCO—A+^li^l + f). (53) 

Since nowadays H and po have small values, we need to choose G(T) = G-(T) and then, using the formula (|4"5j) . 
we get 



m «* (f f j + f (1 - ^1 + f j - a. (54) 

what shows that the effective formulation of LQC is a teleparallel theory. 

To compare with FRW cosmology, since nowadays T is small as compared with p ci we can expand G(T) up to 
second order in T, to get 

G(T) = -A-- + —, (55) 
2 4jO c 
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and inserting this expression into (|45j) . one obtains the first order correction to the FRW cosmology. Moreover, in 
order to obtain exact formulas, one has to use Eq. (14"?} , because in that case G~ l {p) has a very simple, quadratic 
expression, namely 



G-\p) = -2{p + K) (l-^). 



A straightforward calculation yields 



1 , f(p)( Pc -2(p + A)) 3 f(p)(p c -2(p + A)) 

Ueff(P) = -1 + 7 — TT7 ? — TTT, qdec{p) = -1 



(p + A)(p c -(p + A))' 2(p + A)(p c -0» + A))' 

;/ v _ 9 f(p)[-(f'(p) + l)(Pc - 2{p + A)) + 2/(p)] 
JW 2 (p + A)(pc-(p + A)) 

{P> (p + A)(p c -(p + A)) + f(p)(p c -2(p + A)y [ ] 

On the other hand, as we already saw in Sect. II, in brane cosmology in the Randall-Sundrum scenario the modified 
Friedmann equation depicts the following hyperbola 

n P + A (, , P + A 



on the plane (H, p) . 

Finally, by comparing this equation with (52) and making the change p c — > — 2A, it follows that one can view loop 
quantum cosmology as brane cosmology with a negative brane tension and a time-like extra dimension. Therefore, 
in order to obtain the corresponding state-finder parameters formulas in brane cosmology, we just need to do the 
replacement p c — > — 2A. Moreover, with this replacement we can apply the general formulas (|56[) to any EoS, in 
particular to the model studied in Sect. IV. 

A very important remark is here in order. As a result of the above replacement it follows that the dynamics 
resulting for both theories, LQC and BC, are very different, because their corresponding Freedmann equations depict 
two completely different curves. In particular, Rip singularities, as we have seen in Sect. V, are allowed in BC because 
the hyperbola is an unbounded curve. But, since in LQC the Friedmann equation depicts a bounded curve (an ellipse 



Rip singularities cannot appear in this case. For example, for the EoS p(p) = up the universe is non-singular (see 54 1 



for a detailed explanation). In fact, for u > —1 (resp. u < —1) it moves in anti-clockwise (resp. clockwise) way from 
the anti de Sitter solution H = — »J~A^Jl — A (resp. de Sitter solution H = y^t/l — A) towards the de Sitter one 



H = yjl - A (resp. anti de Sitter one H = -y/$Jl - A). 

In the early universe the parameter p c is very large as compared with po (pc 3> po). In principle, p c is fixed 
and has the value p c = 3/(7g J A|g C ), being j BI = 0.2375 the Barbero-Immirzi parameter and X\q C — V^Ibi/^ 
(see, for instance, 52]). But one could also assume that p c can vary with time or, in other words, we can consider 
LQC with the parameter p c in some cosmological epoch. For illustration we can estimate the possible values of p/p c 
at the present epoch (at redshifts < z < 1.75). As the LQC equations formally coincide with the cosmological 
equations on the brane in the RS scenario, one can confront the simplest loop quantum cosmology model (po = A) 
with observational data, in the same manner as it was done in Sect. IV. The analysis shows (see Fig. 6) that the best 
fit for the SNe+i? (z)+BAO data is achieved for 7 = (p + A)/p c m 0.03. Maybe this result can be considered as an 
argument in favor of loop quantum cosmology. 



VII. CONCLUSIONS 



We have confronted in this paper several DE models on the brane with combined data coming from different, 
totally independent cosmological surveys. The analysis here performed shows that the fit of these observational data 
is actually better for the FRW cosmology in the framework of the chosen model for dark energy. Also, owing to 
the fact that the LQC equations formally coincide with the cosmological equations on the brane in the RS scenario, 
we could additionally confront the simplest (but very important) loop quantum cosmology model (po = A) with 
observational data, in the same manner as it was done with the brane models. The analysis we have carried out has 
shown (Fig. 6) that the best fit for SNe+-ff (z)+BAO data is achieved for 7 = (p + A)/p c sa 0.03, what could be 
actually viewed as a good argument in favor of loop quantum cosmology. 
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Figure 6: The 68.3% (solid) and 95.4% (dotted) confidence level contours in the Qa - 7 parameter space from the 
SNe+//(z)+BAO data analysis in the case when po = A. The best-fit parameters for the observational data are Qa = 0.712, 
p/p c = 0.03 with x 2 = 560.58 (to be compared with xLn = 561.31 for the ACDM model). 

Taking everything into account, the observational cosmological results do not exclude, in principle, that the real 
cosmology of the universe we live in could in fact differ from that of the standard FRW model. A window remains 
still open for discrepancy. The importance of joint analysis of the various, independent classes of observational data 
sources has been clearly manifested in the discussion of the different tables and plots. Taking into account together 
SNe apparent magnitude measurements, Hubble parameter evolution data, and BAO and matter density perturbation 
data, we are able to get a quite rigid constraint on the allowed value of the brane tension, in the framework of the 
different brane models considered here. 
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